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Ph Abstract: We derive quantum kinetic equations for scalar fields undergoing coherent evo- 

lution either in time (coherent particle production) or in space (quantum reflection). Our 
central finding is that in systems with certain space-time symmetries, quantum coherence 
manifests itself in the form of new spectral solutions for the dynamical 2-point correlation 
function. This spectral structure leads to a consistent approximation for dynamical equa- 
tions that describe coherent evolution in presence of decohering collisions. We illustrate 
the method by solving the bosonic Klein problem and the bound states for the nonrel- 
ativistic square well potential. We then compare our spectral phase space definition of 
particle number to other definitions in the nonequilibrium field theory. Finally we will 
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explicitly compute the effects of interactions to coherent particle production in the case of 

j-jl an unstable field coupled to an oscillating background. 
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1. Introduction 

Quantum transport effects are gaining more and more interest in many applications in 
modern particle physics and cosmology. This is true in particular for the case of the 
electroweak baryogenesis [1, 2, 3, 4, 5], but also for the leptogenesis [6] or particle creation 
in the early universe [7] and during phase transitions [8] . We have recently developed new 
quantum transport equations for fermionic systems including nonlocal coherence, either 
in space (quantum reflection) or in time (coherent pair production) [9, 10], including the 
effects of decohering collisions [10]. Here we will introduce a similar formalism for the 
scalar fields. As in the fermionic case, the coherence information is found to be encoded in 
new spectral solutions in the phase space of the dynamical 2-point function. The physical 
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information of particle numbers or fluxes and of coherence is carried by a set of scalar 
functions that parametrize the spectral shells in the full 2-point correlation function. 

Our approach can be summarized as follows. We first formulate Schwinger-Dyson 
equations for the 2-point correlation functions using the Closed Time Path (CTP) method. 
We solve the resulting Kadanoff-Baym (KB) equations for the 2-point Wightman func- 
tion in the noninteracting mean field limit in the mixed representation. The most 
general solution for in this limit is a sum of singular spectral distribution functions 
corresponding to the usual mass-shell solutions with a dispersion relation uj'^ = k'^ + m?, 
and a new coherence solutions living at shell /c^ = in the static planar symmetric case 
or at shell A;o = in the case of a spatially homogeneous system. We then turn back to 
the full KB-equations including collision terms and integrate them with an appropriate set 
of moments. On the adiabatic boundaries the lowest moments of A*^ can be related to 
the spectral on-shell functions. From practical point of view the most important aspect 
of our method is that the singular shell structure reduces integrated equations of motion, 
including the collision terms, to a closed set of equations for spectral on-shell functions, 
or equivalently for a finite number of lowest moments of A^ (three in case of a single real 
scalar field). 

Since the singular shell representation of the coherence is so crucial for our formalism, 
we have given several examples which illustrate their physical role. Firstly, we will solve the 
bosonic Klein problem. We will show that in the absence of the coherence shell the quantum 
nature of reflection is completely lost, but that the correct reflection and tunneling factors 
are recovered when the coherence shell is included. We will also show that the spectral phase 
space definition of particle number in our formalism is consistent with other definitions for 
nonequilibrium systems in the literature [11]. In particular, our particle number, when 
applied to Bunch-Davies vacuum in the inflation, corresponds to the adiabatic particle 
number that remains always zero in a conformally coupled scalar theory [12]. The Klein 
problem example also allows us to demonstrate how the on-shell functions are related 
to moment functions that must be used to formulate a dynamical problem with only an 
incomplete information about the variables defining the system. We will also consider 
production of unstable scalar particles by a coherent time dependent background potential 
and decoherence and thermalization of an initially highly correlated state. 

This paper is organized as follows. In section 2 we briefly review the basic CTP- 
formulation for the calculation of the 2-point function and in section 3 we derive the 
spectral shell structure of the Wightman function in the mean field limit. In section 4 
we use our formalism to solve the bosonic Klein problem. We also derive an expression 
for our on-shell particle number in terms of moment functions and compare it with other 
definitions in the literature and apply it to the particle production during inflation. We also 
compute other measurable quantities such as energy density and the pressure. In section 5 
we solve the nonrelativistic problem with a Schrodinger equation and show that there one 
obtains similar coherence solutions for the description of reflection in the planar symmetric 
case. We complete this section by solving the bound states of the square well potential with 
our formalism. In section 6 we show that, just as with fermions, the coherence solutions 
are excluded from the spectral function by the spectral sum rule. In section 7 we derive 
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Figure 1: Schwingcr-Keldysh path in complex time. 

the dynamical (moment) equations for a scalar field including collisions for a spatially 
homogeneous system. In section 8 we consider coherent production of unstable particles. 
Finally section 9 contains our discussion and outlook. 

2. CTP-formalism for scalars 

The main object of interest for us is the 2-point Wightnian function for a real scalar 
field, defined as: 

iA<{u, v) = {cPiv)cP{u)) = Tv{p <P{v)(l){u)} , (2.1) 

where p is some unknown quantum density operator that gives the complete information 
of the system. Instead of trying to find a solution for p, we will set up equations for 
the "in-in" correlation function (2.1) using the Schwinger-Keldysh or Closed Time Path 
(CTP) formalism [13, 14]. In this formalism one first defines 2-point correlation function 
on a complex time-path: 

iAc{u,v) = {Tc[<P{u)(l>{v)]), (2.2) 

where Tc defines time ordering along the Keldysh contour shown in figure 1. One can show, 
for example by use of the two-particle-irreducible (2PI) effective action techniques [14, 15], 
that Ac{x,y) obeys the contour Schwinger-Dyson equation: 

Aciu,v) = A^ciu,v)+ [ dSi / d^Z2 A°ciu,zi)ncizi,Z2)Aciz2,v) , (2.3) 



where A^ is the free propagator of the theory, and the precise form of the self-energy 
function Uq depends on the Lagrangian and the truncation scheme. Once the theory is 
specified, it can be computed from the 2PI-effective action by functional differentiation: 

Uc{u,v)^i /^f^\ , (2.4) 
dAc{v,u) 

where r2 is the sum of all two particle irreducible vacuum graphs in the theory. The 
complex time Green's function in (2.2) can be decomposed in four different 2-point functions 
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with respect to usual real time variable: 

iA<{u,v) =iA+~(n,t;) = {4>{v)<p{u)) 
iA>{u,v) = iA'^{u,v) = {(f>{u)(p{v)) 

iAp{u, v) = iA^^(u, v) = 6{uq — VQ)iA^{u, v) + 9{vq — uo)^A^(u, v) 

iAp{u, v) = iA {u, v) = 9{vq — no)iA^(u, v) + 6{uo — VQ)iA^{u, v) , (2.5) 

where Ap and A^ are the chronological (Feynman) and anti-chronological (anti-Feynman) 
Green's functions, respectively, and A^ and A'* are the Wightman functions we are pri- 
marily interested in solving here. A similar decomposition can be done for the contour 
self-energy Uq to get: 

W'^u, v) = m6-|^4^ ' (2-6) 

where the indices a, b refer to the position of the arguments u and v, respectively, on the 
complex Keldysh time path. When a = -|-1(— 1) the time argument in u belongs to the 
upper (lower) branch in figure 1, and we will use the same notation: 11^ = 11"' , etc. for 
the self energy as we did for the propagators (2.5). It can then be shown that the complex- 
time equation (2.3) is equivalent to the following matrix equation with a usual real time 
argument: 

Ag^ A = 0-3 (5-hn(g)a3A, (2.7) 

where 

and (73 is the usual Pauli matrix, and we defined a shorthand notation for the convolution 
integral: 

f®g = j d''zf{u,z)g{z,v). (2.9) 
We have also left out the labels u and v where obvious; for example 5 = 5'^{u — v). 
2.1 KadanofF-Baym equations 

It's appropriate to further define the retarded and advanced propagators (a similar decom- 
position obviously holds for the self energy function 11) : 

l^'{u,v)= e{u° - v^){A> - A<) 

A'^(n, v) = -6{v^ - u°)(A> - A<). (2.10) 

Moreover, the hermiticity properties of the Wightman functions: 

[iA<'>{u,v)]^ = iA<^>{v,u) (2.11) 

imply that [iA'"(t(, u)]^ = — iA"(f,n), which suggests a decomposition into hermitian and 
antihermitian parts: 

Ah ^^iA^ + A'-) 

i(A'^-A^) = ^(A>-A<). (2.12) 
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The antihermitian part A is called the spectral function. Based on (2.10) it is easy to 
show that Ah and A obey the spectral relation: Ah{u,v) = —isgii{u^ — v^)A{u,v). Since 
the self-energies 11 satisfy identities similar to (2.11), we can define the hermitian and 
antihermitian parts of 11^''" as well: 

r^^^{u^-u^)='- (n> - n<) . (2.13) 

Using the definitions (2.12)-(2.13) it is now straightforward to show that eqs. (2.7), when 
written in the component form, become: 

(Ag 1 -Uh)(S)A-T0Ah = O 

{A^^ - nH)^AH + r^A = d (2.14) 

and 

(Ag 1 -Uh) ^ A< -U< ^ Ah = ^ A< - U< ^ A>) , (2.15) 

where Aq ^ is the inverse free propagator. Equations (2.14) are called the pole equations 
and eq. (2.15) the Kadanoff-Baym (KB) equation. The other KB-equation for A> need not 
be considered, since form the definition (2.12) it immediately follows that A^ = A^ — 2iA. 

2.2 Mixed representation 

The final step in our formal derivation is moving to the mixed representation, to separate 
the external and internal degrees of freedom in the correlators through a Wigner transform: 

F{k,x) = j d\e'^-''F{x + r/2,x-r/2), (2.16) 

where x = {u + v) /2 \s the average coordinate, and k is the internal momentum variable 
conjugate to relative coordinate r = v — u. Performing the Wigner transformation to 
eqs. (2.14) and (2.15) we get the pole-equations 



Ao M - e-^^{n^,}{^} - e~'^{T}{AH} = (2.17) 

A^^Ah - e-''>{IiH}{AH} + e-^'>{T}{A} = 1 (2.18) 

and the KB-equation for A^ becomes 

Aq ^A< - e-'^{nH}{A<} - e-*^{n<}{A^,} = Ccoii , (2.19) 

The collision term in eq. (2.19) is given by 

Cecil = -ie-'"^ ({r}{A<} - {n<}{^}) , (2.20) 

and the 0-operator is the following generalization of the Poisson brackets: 

0{/}{5} = \ [dxf ■ dk9 - dkf ■ dxg] . (2.21) 
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Equations (2.17)-(2.19) are the master equations appropriate for all analysis to be per- 
formed in this paper. Explicit forms of Aq and the interactions depend on the model. In 
this paper we consider a theory defined by the Lagrangian 

^ = lid^<i>f - l^""^^ + Ant , (2.22) 

where m = m{x) is possibly spacetime dependent mass and £int is the interaction part 
to be defined later. The inverse free propagator corresponding to eq. (2.22) in the mixed 
representation is 

Ao ^ ^k^-^d^ + ik-d- m^e-i^"^^ , (2.23) 

where the (9™-derivative always acts on the mass term, and the (9^-derivative to the Green's 
function A, Ah or A^ in eqs. (2.17)-(2.19) respectively. 

3. Shell structure 

In analogy to what was found in the case of fermions [9, 10], a reasonable approximation 
scheme can be developed for an interacting theory employing the spectral shell structure of 
the noninteracting theory. So let us first consider noninter acting fields, for which 11"* = 0. 
In this case the KB-equation (2.19) for A^ decouples from the pole-functions Ah and A 
and reduces to Klein-Gordon equation in momentum space: 

(^k^ -\d^ + ik.d- m^e-'^^^-^^)iA< = . (3.1) 

This is still a very complicated equation because it contains derivative operators to ar- 
bitrarily high orders. We shall analyze it in more detail in the mean field limits where 
m^e~2^™'^fe^ — > m?, assuming also particular space-time symmetries: a case with a spatial 
homogeneity and a static problem with a planar symmetry. 

3.1 Spatially homogeneous case 

In the spatially homogeneous case all spatial gradients of the mass m = m{t) and the 
correlator A^ vanish. Breaking the equation (3.1) into real and imaginary parts and 
expanding to zeroth order in the time derivatives acting on the mass then gives: 

{kl - P -m'^it) - ^dfyA<{k,t) = (3.2) 
kodtiA<{k,t) = 0. (3.3) 

At this point it is relevant to make a comparison with the similar problem involving 
fermions [9, 10]. In the fermionic case the noninteracting KB-equations can be divided 
into algebraic constraint equations which define the phase space structure of the theory 
and to dynamic equations containing all time derivatives. Here such a division is not possi- 
ble, as both equations contain derivative terms, and the shell structure is less obvious than 
with fermions. However, if we assume that /cq 7^ then eq. (3.3) requires that dtiA"^ = 
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at all times and so one must also have d^iA^ = 0. Substituting this back to eq. (3.2) does 
lead to an algebraic equation 

(^kl - P - m^(i)) iA<„, = . (3.4) 
This equation has the spectral solution parametrized by t: 

iA<_,{ko, \k\,t) = 27rsgn{ko)fs,^^{\klt)6{kl -P- m\t)) , (3.5) 
where = sgn(/co). This is just the usual mass-shell dispersion relation 

ko = ±a;g(t) = ±v^|feP + m2(t). (3.6) 

Of course this derivation was not exact, and the solution (3.5) satisfies exactly neither 
(3.2) nor (3.3), except for constant m and fs^g- The point is that equations (3.2) and (3.3) 
are actually inconsistent for a nonconstant mass, but the corrections that would bring the 
consistency back are proportional to mass-gradients. The result (3.5) is thus correct to the 
lowest order in mass-gradients. In particular the effect of the second order derivative term 
in (3.2) to the mass-shell structure is beyond the mean field approximation. 

However, if we first set /cq = 0, then equation (3.3) is identically satisfied and we 
cannot constrain the size of the derivative terms as was done above. Instead, eq. (3.2) now 
becomes: 

52A< = -4a;|(t)A<. (3.7) 
For a constant mass the solution for this equation is 

iA<^,t(fco, \k\,t) = 27rAgcos(2u;^t + 6^) 6{ko) , (3.8) 

where Aj^ and 5j: are some real constants and the (5-function is explicitly taking care of the 
restriction to the shell ko = 0. For a generic time- varying mass an analytical solution for 
zA^ might not be available, but we can write the corresponding solution for A;o = in the 
spectral form: 

iA<{ko,\k\,t) = 27rf,{\k\,mko), (3.9) 

where fc{\k\,t) is some real-valued function. This establishes that there exists a new 
solution living at shell ko = 0. However, one can ask if and how this new solution can 
affect the dynamics of the mass-shell functions f±? Indeed, in the constant mass case, 
where fc is given by eq. (3.8) and f± = const., the answer is no, as expected. To get 
these solutions however, we implicitly introduced prior information on ko that allowed a 
reduction to one particular shell at a time. More generally, one might be interested in 
systems where only an imprecise or even no prior information is available on feg. In such 
cases some integration procedure over ko must be introduced to define observable physical 
quantities, and these quantities typically involve contributions from several shells. When 
such integration procedure is imposed on eqs. (3.2)-(3.3) they generally lead to nontrivial 
mixing involving the functions fc- We will return to this procedure in more detail in section 
4. The basic issue however is that the phase space of the free dynamical function in the 
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noninteracting system is singular in the mean field limit, and it contains new spectral shell 
at /cq = such that the most complete solution for a given momentum | A;| is the combination 
of the solutions (3.5) and (3.9): 

iA<{ko, \k\,t) = iA<_,{ko, \k\,t) + iA<{ko, \k\,t) . (3.10) 

The situation is now seen to be qualitatively equivalent to the case with fermions and we 
interpret analogously [9, 10] that the new /cq = 0-solution (3.9) describes the quantum 
coherence between particles and antiparticles. 

3.2 Planar symmetric case 

Another simple geometry that allows analytic solutions is the case with m = m{z) and 
dt^x.yi^"^ = 0, i.e. a static planar symmetric problem in the average coordinates in the 
Wigner transformation^. In this case the mean field limit of the equation (3.1) is: 

(kl-l? -m'^ + ^dl^i^<{k,z) = Q (3.11) 
kzdJA<{k,z) = 0. (3.12) 

The analysis proceeds analogously to the homogeneous case. For k^ ^ 0, eq. (3.12) gives 
dzi^'" = for all z so that SfiA*" = as well, and eq. (3.11) again reduces to the algebraic 
form: 

(^k^ -P - m2(z)) iA<_^ = . (3.13) 

This equation is similar to eq. (3.4), except that in this case energy is conserved, and the 
mean field momentum kz is the quantity that becomes dependent on z. Taking this into 
account we write the spectral solution in the form 

iA<_g(/co, \kn\,kz,z) = 2irsgn{ko)fs^^{ko, \k^^\, z)6{k'^z - ^m(^)) > (3-14) 

where 

kz = ±kn, = ±^kl-kl-m?{z). (3.15) 

This is again the usual mass shell solution. We get a new solution by first setting kz = 0, 
so that eq. (3.12) is identically satisfied and eq. (3.11) becomes: 

{dl+^kl{z))iK<=^. (3.16) 

In constant mass limit one again finds a solution similar to eq. (3.8), only now restricted 
to shell kz = 0: iKf^^^^ = 27rylcos(2A:m^ + ^) '^(^^)) where ^(/co, | /c|| |) and (5(fco)|^iil) are 
some real constants. For an unspecified spatially varying mass term eq. (3.16) has a generic 
spectral solution: 

iA<(A;o, \k^\\,kz,z) = 2tt fdko, \k^^\, z)5{kz) , (3.17) 

^Note that this actually means that the initial direct space correlator depends only on the internal 
time- or X|| -separation: iA^{t,x;t' ,x') = iA^ {t — t' ; xn — x\^; z, z'). Thus all problems for which the wave 
equations have stationary solutions with ip oc e*^* (such as reflection problem) appear in this sense static 
in the Wigner transformed representation. Same applies of course for stationary dependence on x^. 
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where fc{ko,\k^\\, z) is some real-valued function. Finally, combining the solutions (3.14) 
and (3.17) we find the most complete solution in static, planar symmetric case for given 
energy fco and parallel momentum |A:[||: 

«A<(A:o, \k^^\,k:,,z) = iA<_g(/co, |fc|||, fc^, z) + iA<(A;o, \k^^\,kz,z) . (3.18) 

In this case the mass shell solutions describe modes of left (negative direction in z) and right 
(positive direction in z) moving states, and the zero momentum solution their coherence. 
Indeed, here the coherence solution can be pictured as quantifying the possibility that the 
state is simultaneously going both right and left (interference). It is therefore natural to 
see this solution arising at the average momentum kz = for such a mixture. 

4. Applications 

Having found the spectral solutions, we wish to apply our formalism to solve special physical 
problems including dynamical evolution. Parallel to our analysis with fermions [9, 10], we 
need to define nonsingular weighted 2-point functions to replace the singular ones (3.10) 
and (3.18). These functions can be given in the following generic form: 

f d'^k' - - 

pwiko, k,x) = J -^-^ W{ko, k \ k'^, k'; x) iA<{k'o, k', x) , (4.1) 

where the weight function W{kQ, k \ /cg, k'; x) encodes our knowledge about the energy and 
momentum variables of the state. (In fact a number of different weight functions may be 
needed as we will see below.) Above, with the constant mass examples, we already used 
implicitly weight functions encoding precise energy or momentum resolution and below 
we shall give two more examples. In the case of quantum refiection off a potential wall 
(Klein problem), one may have relatively precise information on energy, but only partial, 
spatially dependent information on the momentum. In case of particle production by a 
homogeneous coherent background field the momentum may be assumed to be known, but 
one has no prior information on the energy, which in this case is not conserved. 

4.1 Klein problem 

As our first example we consider the Klein problem for scalars i.e. a scalar particle reflecting 
off a planar symmetric step potential, see figure 2. We discussed the Klein problem in the 
case of fermions in [9], but the bosonic case has some additional special characteristics. 
This problem could of course be solved by use of the Klein-Gordon wave equation, but 
it provides a nice setting to illustrate how to solve a dynamical problem in our methods, 
which allows us to demonstrate the physicality of our coherence shell solutions. In this case 
the mass is constant, but interactions with the background potential can be represented 
by a singular self-energy correction Hh- It is easy to see that the effect of this correction is 
equivalent to replacing the time derivative with a covariant derivative: idy^ — > idy^ — V{vz) 
in the inverse propagator eq. (2.23). The free KB-equation then becomes 

([fco - V{z)fe-^^^^t -kl + ^dl + ikzdz - m2)iA< = , (4.2) 
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where we have accounted for the fact that the t- and X||-derivatives vanish and we have 
taken fen =0. Taking the real and imaginary parts of eq. (4.2) we find: 

[kl - \dl - kl cos(i5ra,^j)iA< = (4.3) 

(k,d, + kl^\n{\d'^dtS)i^< =0, (4.4) 

where km{z) = \J (/cq ~ ^(■z))^ — with V^z) = VQ{—z). To proceed further we define 
the n-th moments of the Wightman function as integrals over kz'- 

/Ak 
-^k'^iA<{ko,ki=0,kz,z). (4.5) 
Itt 

These integrals are convergent for arbitrary n only if iA^ is compactly supported or a 
sufficiently rapidly decreasing distribution. This is not guaranteed in general. However, 
in this paper we require only that the three lowest moments po,i,2, which are related to 
the free two-point correlator functions, are well defined. Of course, in the adiabatic limit 
where eq. (3.18) is valid, the spectral form of the solution guarantees the existence of all 
moments. These functions correspond to the weighted 2-point functions pw eq. (4.1) with 
specific weig hts Wn = {2TT)''^k^5{ko - A;^)(52(/c'||), which are ah explicitly imposing that 
energy and the momentum parallel to wall are conserved quantities. In the fermionic case 
we only needed one weight function (a 2x2-density matrix) [9] . For the present problem we 
shall need three different functions because of the explicit fc^-dependence in eqs. (4.3)-(4.4). 
Indeed, taking the 0th moment of eq. (4.3) and the 0th and 1st moments of eq. (4.4) the 
following closed set of equations for the moment functions is obtained [16] : 

^d^po + k^mPo - P2 = 
dzPi = 

dzP2-\[dzkl)p^ = Q. (4.6) 

The number of independent moments of course matches the number of independent on- 
shell functions in the spectral solution. Using eq. (3.18) we get the following expressions 
for po,i,2 in terms of the on-shell functions f± and /c for /cq > (and a constant mass m): 

PO = l^{f+ + f-)+fc 

Pi = \u+ - /-) 

P2 = ^(/+ + ^), (4.7) 

while the higher moments are trivially related to pi^2 by (n > 1): p2n+i = k'^pi and 
P2n+2 = k'^~^^p2- Equations (4.6) are our master equations for solving the Klein problem. 
It should be noted that we made no approximations to derive them, since all ^-gradients 
vanish from eqs. (4.6) upon integration over kz- Connection formulae (4.7), and the above 
expressions for the higher moments are formally valid only in the mean field limit. However, 
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Figure 2: Reflection from a step-like potential. Arrows describe the directions of momenta of the 
in- and outgoing particles. 



they can be used to set the physical boundary conditions between the on-shell functions 
f±^c and moments po,i,2 asymptotically in the limit z — > ±00. 

It is easy to solve eqs. (4.6) in the separate regions I and II for the Klein problem 
(see figure 2), where the potential term is either zero or a constant. These solutions will 
contain eight unknown constants that can be fixed by the boundary conditions at 2; ^ itcxD 
and the matching conditions at the potential wall z = induced by the moment equations 
(4.6). First note that all spatial gradients vanish everywhere except at the wall at z = 0. 
Two latter eqs. (4.6) then imply that dzPi.2 = 0, i.e. pi^2 are constants in regions I and 
II. From relations (4.7) it then follows that f± are also constants at \z\ ^0. Now consider 
the boundary conditions appropriate for our reflection problem. The fact that there is no 
incoming flux from the left sets f]^ = 0, because f± were found to be constants. This 
condition also sets coherence solution to zero asymptotically, i.e. /c — > as z — > —00, since 
there are no asymptotic mixing states; note however that we cannot exclude coherence at 
finite distances from the wall based on the boundary conditions alone. Finally, we can 
normalize the incoming flux from the right to unity /I = 1. After these definitions we 
have to account for two distinct possibilities depending on whether the momentum in the 
region II, fej^, is real or imaginary. 

Let us first assume that kg > V + m, so that = q is real. In this case we can have 
nonzero transmitted flux /l^ 7^ 0. Moreover, from the first eq. (4.6) we find that po is 
oscillatory in both regions I and II (fej^^ = k is always real) 

pI, = A' cos(2Cz + 6') + jT^P^ , (4.8) 

where and (5' are new integration constants and P2 are constants related to /j_ by 
eq. (4.7), and we have denoted i =1,11. Combining eqs. (4.7) with eqs. (4.6) we find that 
fc = ~jhrd1pQ, so that coherence solutions are also oscillatory. Since coherence should 
vanish when z — > — 00 we find that A", = in this case. 
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The remaining integration constants are fixed by the matching conditions at z = 
induced by the moment equations (4.6). As dzPi vanishes everywhere we see that pi must 
be continuous over the barrier. Equating f\ = p^^ gives the flux conservation equation: 

/- = 1 - . (4.9) 

The last eq. (4.6) implies that p2 must have a finite discontinuity over the barrier, which 
can be computed by integrating it over a step from z = — etoz = eto give p\ — p^ = 
^[k"^ — q'^]po{z = 0). Using this in the first eq. (4.6) we see that also d^po has at most a 
finite discontinuity over the barrier implying that po and its derivative dzpo are continuous 
at z = 0. Using all these conditions we can fix the remaining integration constants to 
eventually find the transmitted and refiected fluxes: 

(fc - Ai ^ 4fcg 

(fc + g)2' (k + qy 

which are in accordance with the usual Klein-Gordon approach. Finally, the coherence 
solution in the region I can be written in the form: 

fl = l^/ficos{2kz). (4.11) 

It should be noted that if the coherence solution were neglected, the only consistent solution 
for the reflection problem would have been = 0, i.e. that of a classical, complete 
transmission. 

Now consider the case V — m<ko<V + m for which q is imaginary. In this case we 
cannot have mass-shell solutions in region II, so that both f}^ = 0. However, we cannot 
exclude a coherence solution fc = —jhrdlpQ as long as it becomes asymptotically zero as 

m 

z — > — oo. This indeed turns out to be the case: from eq. (4.6) we find that p^^g = and 

p^^ = A^e^l"!^ . (4.12) 

In region I the solutions are of the same form as above above with a real q. We perform 
the same matching procedure over the barrier as in the case of real q to fix the values 
of remaining integration constants. Going through the algebra finally gives the expected 
result with a complete refiection: f\ = 1. The coherence function in the region I is again 
oscillatory: 

fl = -r cos(2kz + 6) with 6 = arcsin ( — ^ — ) , (4-13) 

k -|- \q\ J 

while in the region II it is a dying exponential 

f" = — e^l-?!^ (4 14) 

This vanishes as z — > — oo as required, but remains nonzero in the vicinity of the wall, 
where it clearly describes the quantum tunnelling. Since /j? = 0, the moment function po is 
completely saturated by the coherence fimction. So the tunnelling effect is a pure coherence 
phenomenon that can be interpreted as a maximally coherent virtual pair consisting of a 
left-moving state and its right- moving "antistate". 
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4.2 Particle number 

As another example, we will consider particle number in a spatially homogeneous sys- 
tem. By taking the real and imaginary parts of equation (3.1) subjected to this particular 
symmetry we get now: 

[k'-\d^-m'cos{ldrdi:^))iA< = (4.15) 
(^kodt + m2sin(iar9^j)iA< = . (4.16) 
Analogously to eq. (4.5) we again define the n-th moments of the Wightman function: 

Pnm,t) = I ^k^tA<iko,\klt). (4.17) 
Again, three lowest moments form a closed set of equations: 

^d^ po + uj^po - P2 = 
dtpi = 

dtP2 - ldt{m'')po = . (4.18) 

Note that these equations are exact to all orders of gradients of the mass m{t), assuming 
that the surface terms in k^ vanish. Using the full spectral solution eq. (3.10) one finds the 
following relations of moments to f±^c- 

PO = 7T^(/+ - /-) + fc 

PI = liu + f-) 

P2 = ^(/+-/_). (4.19) 

Unlike the evolution equations, these relations are only valid in the mean field limit. More- 
over, in this approximation, the higher moments are again related to pi^2 by {n > 1): 
P2n+i = "^I^Pi p2n+2 = '^|"''''^P2- Relations (4.19) can be inverted to give f±^c in 
terms of pi, i = 0, 1,2. Following the Feynman-Stiickelberg interpretation we now define 
the phase space number densities for particles and for antiparticles respectively as 

n^ = U{\k\) and = -1 - f_{\k\) . (4.20) 

In terms of the three independent momentum components pi, i = 0, 1,2 we then find: 

1 

= —P2 + Pi 
1 

= —P2 -Pi -I 
^k 

fM) =P0-\p2. (4.21) 
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Using the free theory equations of motion (4.18) we can solve: p2 = + ^'^Po to get 

the expression for the particle number in the form 



= w^Po + ^^tVo + Pi (4.22) 

k 

and the coherence solution now becomes just /c(|^|) = ~c^|/0o/(4ci;|). Note that the moment 
pi remains a constant in a free theory, where dtpi = by eq. (4.18). Setting a constraint 
pi = —1/2 then fixes = at all times, consistent with the fact that (/> is a real scalar 
field. In the operator formalism this constraint is imposed by the Wronskian normalization 
of the mode functions [7]. These results extend trivially to the case of a complex scalar 
field; the only difference is that then can differ from and pi becomes a free parameter 
related to the chemical potential. 

Let us now compare our particle number (4.22) with other definitions in the literature. 
Taking into account the free theory equation of motion (4.18) (and the thermal form of 
the spectral function, eq. (6.6) below), the definition of the particle number by Aarts and 
Berges in ref. [11] can be expressed as 

^k + J)' = PO f^Po + \dlpo] . (4.23) 



k ' 2/ ~ ^"y k^" ' 2 

This agrees with our result (4.22) in the adiabatic, or small coherence limit: d^po « '^gPOi 
as can be directly seen by solving and expanding the square root in eq. (4.23). 

One interesting application is to consider particle number evolution during inflation. 
When applied to expanding space-times in conformal coordinates all that changes in previ- 
ous equations is replacing time with a conformal time: t — > and the mass by an effective 
mass — > = a^m^ — d'^a/a, where a(r/) = —l/{Hri), 77 < 0, is the scale factor. Match- 
ing with the Bunch-Davies vacuimi at early times, one finds that during pure De Sitter 
phase 

PO = -].v\Hl'H-kr,r -° ;^fl + 77^), (4.24) 



2k \ {kr]Y, 

where is the Hankel function of the first kind with u"^ = 9/4 — {m/ H)'^ and H is the 
Hubble expansion rate. It is easy to see that the gradient expansion in the De Sitter case 
can be rewritten as an expansion in 1/|/e77|, so that our inversion formulae (4.19) provide 
a good approximation at early times, but break at the horizon crossing at kr] k, 1. Using 
eq. (4.22), still with pi = —1/2 we find that at early times, or ultraviolet limit k/a ^ 
our particle number behaves as 

This result differs from the particle number defined in ref. [7]: 

1 „9 1 Id/ dr,a 



n 



which at early times times becomes ~ a^^H /2k)'^ . This is not really surprising, because 
the particle number is not unambiguously defined in curved spacetimes. The particle 
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number (4.26) is found by a diagonalization of the Hamiltonian and it corresponds to a 
maximum particle number seen by an ideal detector. Our definition relies on phase space 
arguments and corresponds to an adiabatic particle number [12], which rather tries to 
minimize n. In particular for a conformally coupled scalar theory our particle number can 
be shown to remain zero at all times if it was set to zero in the beginning. 

4.3 Energy density and pressure 

Let us next compute the expectation values of energy (Hamiltonian) density and pres- 
sure, which are the 00- and the ii-components of the energy-momentum tensor (y^'^ is the 
standard Minkowskian metric tensor with the signature (-I-, —,—,—)): 

T'^'' = jr^,d'^<P-Cgf^\ (4.27) 

respectively, in the spatially homogeneous case. For the energy density we get using the 
free theory equations of motion (4.15)-(4.16) 

{Hit)) = {T^\t)) = iUdt^f + Ud^^f + Wcf") 



2 * 

d^fc 



klil^<{ko,\k\,t) 



(27r)3 2 

and for the pressure we get in the same way 

l4 



(27r)4 
d'^A; 1 

(n^ + ng + l)u;g, (4.28) 



(P(t)) = (r-(t)) = j ^^{h^ + kl-uf)iA<{k,^k\,t) 



d^A; 
(2^ 



Pi, _ , / P \ , 

k k 



(4.29) 



There is an explicit contribution from the coherence shell function fc in the pressure, 
signalling that at the quantum level pressure differs from the statistical one. However, 
as discussed in the analysis with fermions [10], we expect that in most cases the direct 
coherence contribution is unobservable in the time-scales longer than the typical oscillatory 
time Atosc ^ 1/uj the coherence solution. However if a strong amplification mechanism 
is in place as in during inflation, even these coherent small scale oscillations might have 
physical consequences. 

Let us point out one delicate issue in computing the energy density and pressure from 
the 2-point function. For example our result for the energy density (4.28) follows directly 
from the direct space integral expression for the Hamiltonian. However, one might instead 
try to start from a partially integrated form of the Hamiltonian density: 

n = l(j) {-df -dl + rP) (j) (4.30) 

Normally this Hamiltonian would give same result as since the only effect of using 
the form (4.30) in (4.28) would be replacing k^ by uj'i in the expression in the second 
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line. Here this difference matters liowever, because in tlie former case the /cQ-sheh does 
not contribute to the energy density whereas in the latter it does. This observation nicely 
underlines the nonlocal nature of our coherence solutions. Indeed, although our equations 
are parametrized by a well defined external time variable, the shell ko = corresponds to a 
completely delocalized, constant mode in the internal variable u — voi the 2-point function 
{(j){v)(j){u)) . Thus, the partial integration and the subsequent neglecting of the boundary 
term leading to the alternative Hamiltonian Ti is not a legitimate operation here. 

5. Non-relativistic case, Schrodinger equation 

Our methods can also be applied to non-relativistic problems. The extension is very simple, 
and we give it here for completeness for a field possibly interacting with some background 
potential V . That is, assume that tp obeys the Schrodinger equation 

idti/j = Hij = ( — —V^ + v]ij. (5.1) 
\ 2m J 

The Wightman function iA^{u,v) = {ijj^ {v)ip{u)) then obeys the equation: 

{id^o + ^Vl-V{u)yA<{u,v)=0, (5.2) 
which in the mixed representation reads: 

s!^^^ + ^L^-^- I'M'-'- = . (5.3) 

This equation resembles the free dynamical equation (3.1) for a relativistic scalar field. The 
physical content of eq. (5.3) can again be most easily analyzed by studying the spatially 
homogeneous and the static, planar symmetric cases. 

In the homogeneous case we assume that the potential and the solutions can only 
depend on time, so that: 

(fco-|^-ni)cos(i5,X))iA<(fc,i) =0 (5.4) 
Q-dt + V{t)sH\dYdtS)i^'^{Kt) = 0. (5.5) 

This set of equations clearly has only positive energy particle solutions with the non- 
relativistic dispersion relation /cq = ^ -|- in the mean field limit, but no negative energy 
antiparticle solutions, nor any coherence solution living at k^ = 0. This was to be expected 
because antiparticles are not automatically a part of the spectrum of a non-relativistic 
field theory. The absence of /cq = 0-shell here is thus consistent with its interpretation as 
describing particle-antiparticle coherence. In static, planar symmetric case one still finds 
the kz = 0-shell solution, which describes the spatial reflection coherence in accordance to 
relativistic fields. 
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5.1 Planar symmetric problems and bound states 

Static, planar symmetric case is more interesting in the nonrelativistic limit. Here the 
equations become identical to eqs. (4.3)-(4.4) for the relativistic scalar field in section 3.2, 
apart from mass-shell dispersion relation, which here reads: 



km = ±\j2m{kQ -V)-kl. (5.6) 

This is of course just the nonrelativistic limit of the dispersion relation used in (4.3)-(4.4). 
Moreover, the moment equations are also obviously identical in form to equations (4.6). 

As an example, let us consider the familiar infinite square well potential in 1-dimensional 
quantum mechanics. Solving the moment equations (4.6) and imposing the same matching 
conditions on the well boundaries as for the Klein problem, one finds that the solution 
consistent with the asymptotic vanishing of po is 



CO „ 

Quiz) 5^ ^n-^ (l + (-1)"+^ cos(2A:„z)) 5{k^ - E„ 



27r 

PO 

n=l 

Pi = 

P2 = eii{z) An—kl 5{ko - En) , (5.7) 

n=l 

where the mass-shell momentum and energy are quantized: km = A;„, = vrn/L and = 
kl/(2m), with n = 1,2, . . .. The function 6'ii(z) = e{z + L/2) - 6{z - L/2) restricts the 
solution to be nonzero only inside the potential well. The remaining constants An can be 
set by the normalization of the solution. We shall soon see that for a pure state these 
constants will correspond to the occupation numbers of the one-particle states labelled by 
quantum number n. Next we shall interpret these results in terms of the spectral shell 
solutions. The generic (mean field) spectral solution for the problem is given by 

iA< = 2n [2mfs,Mkl - kl) + mk,)] . (5.8) 

With this normahzation the relations between fi and the moments read: 

771 

PO = — (/H. + /_)+/e 
Pi = iri{f+- /_) 

P2 = mkm (/+ + /-) . (5.9) 
From equations (5.7)-(5.9) one finds now 

iA< = 9n{z) Y^^nj^ - ^n) + (-1)"""' COs{2knZ)6{k,)\ 2TT6{ko - En) . (5.10) 

n=l 



We can estimate how good approximation this singular shell picture is by computing the 



correlator lAf directly from the one-particle wave functions. For that, we consider a pure 



state 



V / j • 



-17- 






Uj 






-2 









Figure 3: The phase space structure of the correlator iA^ ^{k, z). Shown is the expression inside 
the square brackets in eq. (5.14) as a function of k^/kn at 2 = for ?i = 9 (left) and rt = 99 (right). 



where [a„,am] = Sn,m- The Wightman function corresponding to this state is 

iAl^^{C\i^\t,zmt',z')\0, 
where il;{t, z) = e~^^"'^il)n{z)an is the field operator and 



1pn{z) 



y^2/ L cos{knz)6ii{z) , n = 1,3,5,... 
y^sm{knz)en{z), n = 2,4,6,... 



(5.12) 



(5.13) 



are the normalized one-particle wave-functions in the well. A direct computation of the 
zeroth moment then gives exactly the same expression as in eq. (5.9) with An = fn- 
However, the full Wigner transformed correlator (5.12) becomes now: 



iA<^^{k,z) 



On[z)Y.fJ^ E 

n=l 



± 



sin((/c^=Ffcn)(^-2|z|)) 



+ 



(-1) 



k 



6{ko - En) . (5.14) 



We can see that the entire /c^-dependence of this correlator is encoded into a function 
sin(g(L — 2\z\))/q, which is a representation of the Dirac delta function in the limit kn{L — 
2\z\) » 1: 

sm{q{L - 2\z\)) / q — > Tidiq) . (5.15) 

In this limit the correlator (5.14) reduces to the spectral form shown in eq. (5.10). In 
figure 3 we have plotted the correlator (the expression in the square brackets) eq. (5.14) 
as a function of k^ at the centre of the square well at z = for n = 9 and n = 99 
corresponding to kn{L — 2\z\) « 30 and 300, respectively. It is clear that the phase space 
structure approaches the singular form as the momentum scale kn, compared to the distance 
from the walls, increases. Similar conclusions hold also in the case of mass- or potential 
walls. In particular for a step wall the spectral form for the correlator becomes exact when 
the distance of the wall is large in units 1 /k. 

Let us stress that our only use for the singular shell approximation is to relate the 
moment functions to on-shell functions in adiabatic regime and to provide a practical 
scheme to evaluate the collision term that gives rise to a closed set of equations for the 
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moments po,i,2 (or equivalently the on-shell functions f±^c)- The results of this section 
show that that this scheme can be useful even in most extreme situations. Indeed, even a 
correlator shown in the left panel of figure 3 should be reasonably well represented by a 
singular ansatz in the collision integral whenever the matrix element is a relatively smooth 
function of momentum. For smooth wall profiles and for slowly varying driving forces 
the correlator can approach the singular ansatz even inside the transition region. The 
quantitative measure for this is the ratio of the momentum to the rate of change of the 
potential. Finally, let us point out that our results show that the new kz = 0-shell is equally 
well visible and concentrated as the standard mass-shells. That is, our approximation for 
the phase space is in no means less rigorous than the standard derivation of the ordinary 
Boltzmann equations relying only on the mass-shell contributions. 

6. Spectral function and thermal limit 

Having shown that the coherence solutions are part of the spectrum of the dynamical 
2-point function, let us now show that they do not appear in the pole functions, or in 
particular in the spectral function A. We shall only consider the spatially homogeneous 
noninteracting case. With no interactions the equation for A eq. (2.17) becomes identical 
with the K-G equation (3.1) for A^. Consequently, the most general solution fulfilling the 
mean field phase space constraints must be identical to eq. (3.10) with three yet undefined 
on-shell functions and f^. However, the spectral function must also satisfy the spectral 
sum-rule, which in direct space reads 

2 dt'A{t', u; t, v)t'=t = -iS^{u - v) . (6.1) 

This follows from the canonical equal time commutation relations of the scalar field 4>, or it 
can be derived from the pole equations (2.17)-(2.18). Transforming eq. (6.1) to the mixed 
representation gives 

/ ^{ko+'-dt)A{k,x) = l. (6.2) 

The time-derivative appearing in this representation is usually omitted in the literature, 
by an implicit assumption of translational invariance. This is appropriate for example for 
thermal equilibrium systems, but for more general nonequilibrium problems it does give a 
new independent constraint. In terms of moment functions p-^ = J ^ kQ A the sum rule 
(6.2) becomes: 

Pt = l, dtp^ = 0. (6.3) 

The latter constraint implies that Pq- = const, and furthermore dfp^' = 0. Dynamical 
equations for the moment functions are of course identical to eqs. (4.18) for pn- Imposing 
the sum-rule constraints on these equations we find that also p2' = const. One then finds 
that either Pq^ = P^ = 0, or the mass is a constant everywhere. To get a continuous 
constant mass limit we must always set: 

p{f = p^ = 0, pt = \. (6.4) 
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Connection formulae identical to eq. (4.19) between p„ and f^^ then give: 

ft = \, ff = 0. (6.5) 

These will finally reduce A to its standard thermal form 

A = 7rsgn(/co) 6 {k^ - m^) . (6.6) 

Let us finally consider the thermal equilibrium limit for the function fA^. The new con- 
straining element here is the Kubo-Martin-Schwinger (KMS) boundary condition 

A>(t) = A<(t + z/?) . (6.7) 

However, first note that the relation A^ — A^ = —2iA sets 

/±-/±=l and />-/< = 0. (6.8) 

Then the momentum space version of the KMS-condition: A^(A;o) = e^'^oA^(A;o) is enough 
to set the mass-shell distributions to the statistical limit: 

fko = '^oq(A;o) and =1 + ncq{ko) , (6.9) 

where ncq(A;o) = l/{e^^° — 1) is the usual Bose-Einstein distribution. Moreover, coherence 
functions are subjected to constraint 

/< = /> . (6.10) 

However, the KMS-condition only makes sense when the system has a time-independent 
Hamiltonian. This implies time-translational invariance in real time which immediately 
eliminates coherence, leading to the standard thermal expressions [5]: 

zA<q = 27rsgn(A:o)ncq(A:o)(5 (/c^ - m^) , 

iA> = 27rsgn(A:o)(l + n^^{ko))6 [k^ - m^) . (6.11) 
7. The case with colHsions 

We now move to consider the case with collisions. As explained in the introduction we 
are using exact forms of the integrated evolution equations except for the evaluation of 
the collision terms. It is only there that we need to use the (in general approximate) 
connection formulae (4.19) between the moments and the on-shell functions to get the 
collision integrals in closed form. The complete set of Schwinger-Dyson equations of the 
interacting theory (2.17)-(2.19) are too complicated to be used in practical applications 
without approximations. Here we make a series of approximations that will enable us to 
consider the essential quantum dynamics in terms of the three lowest moments po,i,2 in the 
presence of collisions. 

First, we will consider a weakly interacting theory, so that the usual quasiparticle 
approximation applies. This means that the term oc H/^A^ is included to modify the 
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dispersion relations in both the pole equations (2.17)-(2.18) and in the dynamical equation 
(2.19) for the Wightman function iA'^. However the terms oc T, which are the source of 
broadening of the phase space of the pole functions are neglected in the pole equations 
entirely to allow a singular phase space structure, as usual in thermal field theory. Using 
the constraint — A^ = —2iA one can show that within the quasiparticle approximation 
it is consistent to drop the term oc 11^ A// in the dynamical equation (2.19), and neglect 
the collision term when working out the spectral structure for A^ , even if in the coupling 
constant expansion the dropped terms are of same order as Hh and collision integral C. 
A more complete derivation of these approximations and discussion of the role of different 
self-energy functionals could be found in ref. [5]. Second, we will compute the collision 
term in the r.h.s. of eq. (2.19) only up to first order gradients; this should be a good 
approximation at least for cases where the variations are affecting only a small subset of the 
entire interacting system. With these approximations one can find the desired reduction 
to three moments only. However, we shall here neglect also the term ~ Hj^A^, which 
would just change the dispersion relations of the states, without altering the qualitative 
aspects of collisions on the evolution of the system. This leaves us with the flow term of 
the free theory. The final form of the dynamical equation with collisions in the spatially 
homogeneous case then reads: 

- ^df + ikodt - m2e"5^*"^^^o)iA< = iCcoii, (7.1) 

where the collision term is 

Ccou = - (r iA< - iU<A) + iO ({r}{iA<} - {iU<}{A}) 

= Co + iCi . (7.2) 

Taking the real and imaginary parts of eq. (7.1) gives the coupled equations: 

Co . (7.3) 

Our first task is to find the singular shell structure for zA^. As explained above, consistency 
with the quasiparticle limit of the phase space requires neglecting the collision terms and 
also the gradients of the mass (as explained in section 3), so that we are left with the 
same mean field constraint equations (3.2)-(3.3) as in the free theory case. We then find 
the familiar shell structure for the Wightman function: 

iA< = 27r (sgn(A;o)/.,/(A;2 - m^) + fc6{ko)) , (7.4) 

and the same relations between the functions f±^c and the moments po,i,2 as in the free field 
case, given by eq. (4.19). These relations are the core of our approximation scheme, since 
they allow the equations of motion derived from (7.3) to close with only the three lowest 
moments. Indeed, integrating both equations in (7.3) with a flat weight and the second 



(Mt + m2sin(iar<))iA< = 
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equation weighted by ko, we get the following generalizations of the free- field moment 
equations (4.18): 



1 



d?Po + w|/0o - P2 = - (Ci) 



4 




(7.5) 



where the collision integrals appearing on the r.h.s. of eqs. (7.5) are 




(7.6) 



The problem with these equations is that functions T and 11^ can have an arbitrary phase 
space structure, so that the collision integrals are a priori not related to the moments pi in 
any simple way. That is, equations (7.5)-(7.6) do not close. This is of course to be expected, 
because integration erases a lot of information from the system. Equations (7.5)-(7.6) are 
in fact useful only if the collision terms can be reasonably well approximated by some 
expansion in the lowest moments. This is precisely what our singular shell structure for 
the Wightman function iA^ does. Indeed, when the structure (7.4) is fed into the collision 
integrals (7.6), they become completely parametrized by the on-shell functions /±,c, which 
on the other hand are related to the lowest moments po,i,2 via eq. (4.19). Note that this 
approach is more elaborate than a simple truncation of the moment expansion, because 
the singular shell structure provides nontrivial information about the phase space of the 
collision integrals. This is of particular importance for the coherence shells, as we shall see 
below. 

To be specific, let us assume a simple thermal interaction for which the self-energies 
do not depend on and obey the KMS-relation 11^ = e^^'^Il^. Moreover, it is natural to 
require (at least in the vicinity of the mass-shell) that T{—ko) = —T{kQ) and T{kQ = 0) = 0. 
These assumptions should hold quite generally for a thermal F; in the appendix A we will 
show explicitly that they hold in the case of a three body Yukawa interaction. Then, 
using the relation lAf^ = 2neqA given by eqs. (6.6) and (6.11), and the inverse relations of 
eq. (4.19) we find: 



(Ci) 




7^ — dko^m{f+ - f~) + dkoTo fc 



'k 



(Co) 



1 



r„[(/+ + /_)-(/:'' + /!^)] 



k 




(7.7) 
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where we have defined = neq(±u;^), and we have neglected terms of order 0(T , dt{m )T). 
The Fj-functions involve projections onto the mass- and the zero-momentum shells: 

Tm{\k\,t) = T{ko = UJ^{t),\k\) 

dk^Tom = dkor{ko = 0,\k\). (7.8) 

Note in particular that the derivative term is computed at "ofi-shell" value A;o = corre- 
sponding to the coherence shell. Expressing /-functions in terms of the moment functions 
through eqs. (4.19), and inserting the resulting expressions back to equations (7.5) we 
eventually find a closed set of equations: 

^d^Po + ^Ipo - P2 = -^dk^To dtpo 

dtPl = -—Tm {Pl - Pl,cq) 
dtP2 - -dt{w?)po = r„, {p2 - P2,cq) • (7.9) 

^ ^k 

eqs. (7.9) are the master equations that are used in section 8 to study the coherent pro- 
duction of unstable particles in an oscillatory background. 



8. Coherent production of unstable particles 

In this section we shall compute the effect collisions on the coherent production of unstable 
scalar particles in the presence of a time varying driving field. This problem is very similar 
to the one we considered for fermions in ref. [10]. In fact we are taking the mass term 
driving the particle production to be of the same form as in [10]: 

m^{t) = \mo + e-'^^{Acos{2uJ^t) + iB sm{2uj^t))\^ , (8.1) 

where mo, A, B, uj^p (oscillation frequency of the the driving field ip) and 7 are real con- 
stants. To illustrate more clearly the qualitative aspects of the method we take some 
parameters of the model in this example to be outside the adiabatic limit. Reader is 
warned that this may make quantitative results somewhat inaccurate. It should be em- 
phasised that the method, especially the calculation of is proven only in adiabatic limit. 
The task is simply to compute the collision functions Tm and dk^To appearing in (7.8) for 
the particular model under consideration, and solve the equations (7.9) with some suitable 
initial conditions. Here we shall consider the interaction 

Ant = -y V^V"/-, (8.2) 

where (j) is the real scalar field whose dynamics we are interested in, and is some fermion 
field, which will be assumed to be in thermal equilibrium at all times. We assume that 
m0 > 2m^ at least for some t, and consider the effect of the induced instability on the 
(/)-particle production. The explicit expressions for Tm and dk^TQ that we will be using are 
computed in the appendix A. 
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Numerical solution of equations (7.9) is not always straightforward; depending on the 
precise form of the driving term they can become very unstable against numerical errors. 
The problem can be traced to the third equation in (7.3). It turns out that for strong 
driving fields the dynamical evolution of p2 is very delicate and it is impossible to discern 
the physical solution from exponentially growing spurious numerical errors. However, these 
instabilities can be circumvented by transforming equations (7.9) into an equivalent set of 
nonlinear equations in the presence of a constraint. Indeed, taking the sum of the first 
equation in (7.9) multiplied by —2dtpo and the third equation multiplied by 2pQ we will 
get first order differential equation: 

dtX = -2dtpo {iCi) + 2po {koCo) , (8.3) 

where 

X = 2p2Po-u;lpl-\{dtPo)'. (8.4) 

One observes that the function X is a constant of motion in free theory, and so it should 
only change slowly in the interacting theory providing collision integrals are small, which 
should be the case for the perturbation expansion to be valid. Initial value of X can of 
course be calculated from initial values of pi'.s. The advantage of this formulation is that 
we can use the algebraic equation to solve p2, while replacing the dynamical equation for 
P2 by a much better behaving equation for X. That is, we use the equations 

^dtY + u^po - p2 = -]^dkoTo Y 

dtPl = (pi - pi cq) 

dtX = dk^To y2 _ 2^^E!!l - p2,cq) , (8.5) 

^k 

where Y = dtpo and p2 = {X + + uip^) /2pQ. Moreover, thermal values for pi^2 can be 
seen from eqs. (4.19) and (6.9)^: pi^eq = —5 and P2,eq = ^k{''^eq{^j^) + |)- Equations (8.5) 
are numerically stable and easy to solve. In figure 4 we show the evolution of the number 
density and the absolute value of the coherence function fc{\k\) as given by eq. (4.21): 

ng=— p2-^ and fc{\k\) = po-\p2, (8.6) 
^k ^ 

under the influence of a driving mass term (8.1). The uppermost panel corresponds to 
case without collisions. The increase of the number density (thick solid line) is seen to 
be accompanied by a steady growth of the amplitude of the coherence (rapidly oscillating 
thin dotted blue line). In the middle panel we show the same solution in the case where 
we have included only the mass-shell collision terms, but set artificially dk^TQ = 0. Here 
one sees that the number density decreases in the intervals between the stepwise growth 

^The equation for pi in eqs. (8.5) is actually trivial in the case of a real scalar field where pi = — | 
throughout. However, these equations are valid as such also for a complex scalar field, where pi becomes a 
dynamical variable related to a chemical potential. 



- 24 - 




2 4 6 8 10 12 14 16 

X 



Figure 4: Shown is the number density rij: (thick sohd hne) and the coherence function /c(|fc|) 
(dotted(blue) hne) with changing interactions. The driving mass function is taken to be ra^it) = 
1(1 + 1. 5cos(2a;^t) + i0.1sin(2w^i))rp. The upper panel corresponds to the case without coUisions. 
In the central panel we have included the collision terms on the mass-shells, but kept dko^o = 0. 
In the lowest panel the full interaction terms were kept for all shells. For parameters we have used 
|fc| — 0.6 T, y = 1, — 0.1 T and uj^p — O.IT, where temperature T sets the scale. At t = lJJ^pt — 
the system is in the adiabatic vacuum. 



as a result of decays. Finally the lowest panel shows the case where all collision terms are 
included properly. It is interesting that the growth of the number density is most strongly 
affected by the collisions acting on the coherence solution. This can be understood when 
one solves for the time-evolution of directly from (7.9). The result is: 

dtnj: = {dtuj^)fc - ^Jm{n^ - n r) . (8.7) 
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That is, the creation rate of is completely controlled by the coherence. Thus collisions 
that destroy the coherence also directly cut down the growth of the number density. This 
is potentially more important effect than the decay of the already created particle number 
effected by the on-shell decay rate Tm- 

From equations (8.5) it is not obvious how collisions affect the coherence. However, 
it can be shown that their only stationary solution corresponds to l^q = and X^q = 
(neq(w^) + |)^, and this equilibrium limit can only be reached only through the effect of 
SfegTo-terms. That this solution corresponds to vanishing coherence, becomes most evident 
when the evolution equations are written directly in terms of fc- The general form of 
eqs. (8.5) is quite messy when written in terms of the /j-functions, but they become much 
simpler if we assume that the driving mass term is a constant. In this limit the growth 
term for vanishes in eq.(8.7) and the equation for fc becomes: 

52 /c + 4L^fc/c = -2dk,rodtfc . (8.8) 

That is, when the driving term is shut off, the particle number decays towards the equilib- 
rium value with a rate given by Tm as expected. Simultaneously the coherence is driven 
to zero by the interaction term dkoTQ. Note that collisions act as a friction on coherence, 
just as in the case with fermions [10]. Thus the solution is an oscillating function with 
an exponentially dying amplitude. We demonstrate this behavior from the full equations 
in figure 5 for a particular initial configuration created from vacuum by the same driving 
term used in figure 4. 

Finally, in figure 6 we show the evolution of the number density under the influence 
of a driving force whose amplitude decays in time with the exponent 7 = 0.05. The black 
solid line corresponds to a noninteracting case, and the thin solid (green), dash-dotted 
(blue) and dashed (red) lines to interacting cases with increasing strength of interaction. 
The dotted line corresponds to the local equilibrium number density Ueq- This function 
oscillates as a result of the time-varying mass of the field. The steady growth of the 
particle number is now absent, as a result of the non- resonant nature of the driving field. 
The instantaneous particle number is mostly controlled by the quantum coherence effects 
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Figure 6: Number densities of particles in same configuration as before but with |fc| = 0.3 T and 
mass oscillation amplitude suppressed by exponent with 7 = 0.05. The black solid line corresponds 
to a noninteracting case y = 0, and the thin solid (green) y ~ I, dash-dotted (blue) y = 3 and 
dashed (red) y ~ 5 lines cases with different interaction strengths. The thin dotted line corresponds 
to the local equilibrium number density ricq. 

when the interactions are absent or weak. In the case with weakest interactions the particle 
number is still strongly modulated by the background field, but there is a clear decaying 
trend due to interactions. In the case shown by dash-dotted (blue) curve the interactions 
are so strong that they almost eliminate the coherence after the first peak and the tendency 
of interactions to push the particle number towards their equilibrium values is beginning 
to show. This tendency is even clearer in for the dashed (red) curve corresponding to 
most strongly interacting case. In the limit of infinitely strong interactions no coherence 
would be left and the particle number would follow the local adiabatic equilibrium particle 
number. 

These examples show that using our methods it is possible to describe coherent particle 
production in the presence of decohering interactions. We considered only the case of 
decay, but it would be straightforward to extend the treatment to the case of collisions. 
The effect of collisions was qualitatively different on the mass-shell and on the coherence 
solutions, since in the latter case interactions introduce a friction term that tends to erase 
the oscillating coherence function, whereas the mass shell distributions were found to feel 
the usual relaxation towards equilibrium. 

9. Conclusions and outlook 

In this paper we have derived quantum transport equations including nonlocal coherence 
effects for relativistic and nonrelativistic scalar fields in systems with particular spacetime 
symmetries. The time-dependent but homogeneous systems include for example particle 
production during phase transitions or during the inflation in the early universe. The 
static, planar symmetric problems include for example reflection off a potential, such as 
Klein problem, or off a mass wall induced by a phase transition front in the early universe. 
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The key observation leading to a calculable approximation scheme was the observation 
that in both of these geometries the 2-point correlator, written in the mixed representation, 
has new spectral solutions living on shell /cq = in the homogeneous case and on shell 
kz = in the static planar symmetric case. These solutions were interpreted to carry 
information on the nonlocal coherence between particle and antiparticle excitations in the 
former and between left- and right moving states in the latter case. We demonstrated 
the physical nature of these new spectral solutions by applying the formalism to exactly 
solvable models, such as relativistic Klein problem and bound states in one dimensional 
nonrelativistic potential well. 

The nontrivial singularity structures described above were found as an approximation 
to collisionless equations to the lowest order in gradients. The core of our calculational 
scheme was the argument that these structures should provide a reasonable ansatz for the 2- 
point function when relating the moments to physical on-shell functions at the boundary of 
the system and when computing collision terms in the moment expansion of the Kadanoff- 
Baym equation. When the ansatz is introduced into the collision terms they become 
completely parametrized by the on-shell functions in the ansatz. Because the on-shell 
functions can be uniquely related to the lowest moment functions, the resulting moment 
equations can be solved. Despite the simplicity of the resulting equations, they contain the 
essential information about the evolution of nonlocal quantum coherence under decohering 
interactions. Based on the nature of the approximation we argue that our method could be 
useful even when the background field is not necessarily slowly varying. Method requires 
the existence of an adiabatic boundary, or boundaries where the on-shell functions can 
be related to the moment functions. It is thus best suited to problems with localized 
disturbances in background field configurations with asymptotic adiabatic boundaries. 

Our method provides a natural definition for the adiabatic particle number related 
to the value of the phase-space functions multiplying the singular mass-shells. This defini- 
tion was applied to the particle number evolution during inflation where it was shown to 
correspond to the adiabatic particle number defined e.g. in ref. [12]. Moreover, our particle 
number coincided with the definition of ref. [11] for slowly varying fields. We also computed 
the particle number evolution in the presence of a driving background interaction in the 
form of a time-dependent mass term. This situation could model for example the particle 
creation during phase transitions or at the end of the inflation. We then included deco- 
herence assuming that the scalar particles created by the background fields were unstable. 
We found that the effect of interactions on particle number divided to two parts: first the 
existing particle number was suppressed by decays as expected and secondly interactions 
provided a friction term on the growth of the coherence. However, as the growth of was 
found to be completely controlled by coherence, the friction term turned out to be most 
efficient in reducing the particle number created by unit time. 

Many of the results presented here were qualitatively similar to those derived earlier 
by us [9, 10] for fermions. However, details of the derivation were substantially different 
so as to warrant a complete independent treatment in this paper. It would be interesting 
to apply our formalism to study for example the effect of collisions or decays on the pro- 
duction of scalar fields in a realistic model for a parametric resonance. It should also be 
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straightforward to extend our nonrelativistic formulation for example for 3D-cubic lattice 
potentials and study atoms in such lattices under the influence of external thermal noise. 
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A. Yukawa interaction with thermal background 

In this appendix we compute the appropriate self-energies for the use of the master equa- 
tions (7.9) in the case of a Yukawa interaction with thermal background. We start with 
the interaction Lagrangian: 

Ant = -y 4^i^<P (A.i) 

where (f) is the considered real scalar field and "0 is some fermion field. We use the 2PI 
effective action method to calculate the self-energies (2.4) at 1-loop level. The lowest order 
2PI-graph based on interaction (A.I) (figure 7) gives the contribution 

r2Pi = -y2 j d\ d\ Ti [G{u, v)G{v, u)] A{u, v) , (A.2) 

where G is the propagator of the fermion ip and the integration is over the Keldysh time 
path. From this we get the self-energies by use of eq. (2.6). In particular, after performing 
the Wigner transformations we have: 

m<'>{k,x) = -y^ I ^^Tr[G>'<{k',x)G<'>{k + k',x)] . (A.3) 

We assume thermal background so that the fermion ip distributions appearing in the loop 
are thermal. The appropriate thermal propagators with real constant mass are (see for 
example [5]): 

iGfq^{k) = 2^ sgn{ko) {jt + mf)n{.^{ko) 5 {k"^ - mj) 

iG>^{k) = 2^sgn(fco) (H^f) (1 - ril^{h))S{k^ - mj) , (A.4) 

where neq(A;o) = i/{e^^° + 1) is the standard Fermi-Dirac distribution function. 

In our present analysis we need to evaluate the self-energies 11^'^ both on the mass- 
shell k^ — K^ = m{t)'^ as well as on the fco = 0-shell. On the mass-shell we get: 



^^^(A„)y^_^ (e. + l)(eVT-. + i)' 



iIi<{ko = ±uj^^{t),\k\)='j^^9{\^) / dy ^-^^^^^3^^^;^-;^ , (A.5) 



with 

a = ^ and /? = ^M, (a.6) 

where Am ^ A(m^, mj, m^) = m?{m'^ — Amj) is the usual kinematic phase space function 
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Figure 7: The only diagram contributing to the self-energy at 1-loop level. 



on the mass-shell. On the ko = 0-shell we get instead: 

where now Aq = A(— = P(P -|-4mj). 

Since we are computing n"'"s in the thermal limit, the expression for 11^ can be 
obtained from that for 11^ by use of the Kubo-Martin-Schwinger (KMS) relation: 

n>(A;) = e^'=''n<(A;), (A.8) 

This relation follows directly from the corresponding relation between the thermal equilib- 
rium propagators Geq^ eq. (A. 4). Using this relation we find that 

r(A;) = ^(l + e^'=o)n<(A:). (A.9) 

Now a direct computation shows that 

r(-A;o,|^|) = -r(A:o,|^|) and r(fco = 0, |^|) = , (A.IO) 

such that the assumptions made in section 7 are indeed verified for this type of interaction. 
For the use of the master equations (7.9) we still need the derivative dkoTQ{\k\) = dkQT{kQ = 
0, |A;|). A direct computation gives 

dk,T{ko = 0, \k\) = ^m<(fco = 0, \k\) . (A.ll) 

Expressions (A. 5) and (A. 7) together with the relations (A.9)-(A.ll) complete the compu- 
tation of all required self-energy functions for the use of the master equations (7.9). 



References 

[1] A.G. Cohen, D.B. Kaplan, and A.E. Nelson, Progress in electroweak baryogenesis, Ann. Rev. 
Nucl. Part. Sci. 43 (1993) 27 [hep-pli/9302210]; 

V.A. Rubakov and M.E. Shaposhnikov, Electroweak baryon number non-conservation in the 
early universe and in high-energy collisions, Usp. Fiz. Nauk. 166 (1996) 493, [Phys. Usp. 39 
(1996) 461] [hep-ph/9603208]. 



- 30 - 



[2] J.M. Cline, M. Joyce and K. Kainulainen, Supersymmetric electroweak baryogenesis, 
J. High Energy Phys. 07 (2000) 018 [hep-ph/0006119]; 

J.M. Cline, M. Joyce and K. Kainulainen. Supersymmetric electroweak baryogenesis in the 
WKB approximation, Phys. Lett. B 417 (1998) 79, [Erratum ibid. 448 (1999) 321] 
[hep-ph/9708393]; 

J.M. Cline and K. Kainulainen, A new source for electroweak baryogenesis in the MSSM, 
Phys. Rev. Lett. 85 (2000) 5519 [hep-ph/0002272]. 

[3] M. Joyce, T. Prokopec and N. Turok, Nonlocal electroweak baryogenesis. Part 2: The 
Classical regime, Phys. Rev. D 53 (1996) 2958 [hep-ph/9410282]; 

Electroweak baryogenesis from a classical force, Phys. Rev. Lett. 75 (1995) 1695 [Erratum 
ibid. 75 (1995) 3375] [hep-ph/9408339]; 

Nonlocal electroweak baryogenesis. Part 1: Thin wall regime, Phys. Rev. D 53 (1996) 2930 
[hep-ph/9410281]. 

[4] K. Kainulainen, T. Prokopec, M.G. Schmidt and S. Weinstock, First principle derivation of 
semiclassical force for electroweak baryogenesis, J. High Energy Phys. 06 (2001) 031 
[hep-ph/0105295]; 

Semiclassical force for electroweak baryogenesis: Three-dimensional derivation, Phys. Rev. D 
66 (2002) 043502 [hep-ph/0202177]; 

K. Kainulainen, T. Prokopec, M.G. Schmidt and S. Weinstock, Quantum Boltzmann 
equations for electroweak baryogenesis including gauge fields, [hep-ph/0201293]. 

[5] T. Prokopec, M.G. Schmidt and S. Weinstock, Transport equations for chiral fermions to 
order h bar and electroweak baryogenesis. Part 1., Ann. Phys. (NY) 314 (2004) 208 
[hep-ph/0312110]; 

Transport equations for chiral fermions to order h-bar and electroweak baryogenesis. Part LL., 
Ann. Phys. (NY) 314 (2004) 267 [hep-ph/0406140]. 

[6] W. BuchmuUer, R.D. Peccei and T. Yanagida, Leptogenesis as the origin of matter, Ann. 
Rev. Nucl. Part. Sci. 55 (2005) 311 [hep-ph/0502169]; 

W. BuchmuUer, P. Di Bari and M. Plumacher, Leptogenesis for pedestrians, Ann. Phys. (NY) 
315 (2005) 305 [hep-ph/0401240]. 

[7] B. Garbrecht, T. Prokopec and M.G. Schmidt, Particle number in kinetic theory ; 
Eur. Phys. J. C 38 (2004) 135 [hep-th/0211219]. 

[8] J.H. Traschen and R.H. Brandenberger, Particle production during out- of- equilibrium phase 
transitions, Phys. Rev. D 42 (1990) 2491; 

Y. Shtanov, J.H. Traschen and R.H. Brandenberger, Universe reheating after inflation, Phys. 
Rev. D 51 (1995) 5438 [hep-ph/9407247]. 

[9] M. Herranen, K. Kainulainen and P.M. Rahkila, Towards a kinetic theory for fermions with 
quantum coherence, Nucl. Phys. B 810 (2009) 389 [arXiv: 0807. 1415]. 

[10] M. Herranen, K. Kainulainen and P.M. Rahkila, Quantum kinetic theory for fermions in 
temporally varying backgrounds, J. High Energy Phys. 09 (2008) 032 [arXiv : 0807 . 1435]. 

[11] G. Aarts and J. Berges, Nonequilibrium time evolution of the spectral function in quantum 
field theory, Phys. Rev. D 64 (2001) 105010 [hep-ph/0 103049]. 

[12] N.D. Birrell and P.C.W. Davies, Quantum fields in curved space, paperback edition, 
Cambridge University Press, New York USA, (1994). 



- 31 - 



[13] J.S. Schwinger, Brownian motion of a quantum oscillator, J. Math. Phys. 2 (1961) 407; 
L.V. Kcldysh, Diagram technique for nonequilibrium processes, Zh. Eksp. Teor. Fiz. 47 
(1964) 1515 [Sov. Phys. JETP 20 (1965) 1018]. 

[14] K.C. Chou, Z.B. Su, B.L. Hao and L. Yu, Equilibrium and Nonequilibrium Formalisms Made 
Unified, Phys. Kept. 118 (1985) 1. 

[15] J.M. Luttingcr and J.C. Ward. Ground state energy of a many fermion system. 2., Phys. 
Rev. 118 (1960) 1417; 

J.M. Cornwall, R. Jackiw and E. Tomboulis, Effective Action for Composite Operators, Phys. 
Rev. D 10 (1974) 2428; 

E. Calzetta and B.L. Hu, Nonequilibrium Quantum Fields: Closed Time Path Effective 
Action, Wigner Function and Boltzmann Equation, Phys. Rev. D 37 (1988) 2878. 

[16] P.F. Zhuang and U.W. Heinz, Equal-time hierarchies for quantum transport theory, Phys. 
Rev. D 57 (1998) 6525 [hep-ph/9801221]. 



- 32 - 



